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Abstract. Jacobi polynomials are mapped onto the continuous Harm polynomials by the 
Fourier transform and the orthogonality relations for the continuous Hahn polynomials then 
follow from the orthogonality relations for the Jacobi polynomials and the Parseval formula. 
In a special case this relation dates back to work by Bateman in 1933 and we follow a part 
of the historical development for these polynomials. Some applications of this relation are 
given. 



1. Introduction and history 

In Askey's scheme of hypergeometric orthogonal polynomials we find the Jacobi polyno- 
mials and the continuous Hahn polynomials; see Askey and Wilson [5, Appendix] with the 
correction in [2], Koekoek and Swarttouw [20] or Koornwinder [23, §5] for information on 
Askey's scheme. In the hierarchy of Askey's scheme of hypergeometric orthogonal polyno- 
mials the continuous Hahn polynomials are above the Jacobi polynomials since they have 
one extra degree of freedom. In this paper we consider a way of going up in the Askey 
scheme from the Jacobi polynomials to the continuous Hahn polynomials by use of the 
Fourier transform. This method is a simple extension of some special cases introduced by 
Bateman in the 1930's. 

In his 1933 paper [9] Bateman introduces the polynomial F n satisfying 

(1.1) ^ n ^\d~ ) cos ^ 1_1 x = cos h _1 £-Pn(tanhx), 

where P n (x) = Pj^ ,0 \x) is the Legendre polynomial, cf. (2.3) for its definition. This 
is possible since 4- cosh - = —cosh - tanhx and ^-tanhx = 1 — tanh x. Bateman 

1 ax ax 

derives in [9] various properties for the polynomials F n , such as generating functions, 
explicit expressions as hypergeometric series, three-term recurrence relation, difference 
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equations, integral representations and more, see also [10]. One year later Bateman proves 
the orthogonality relations, cf. [11], [10], 



(1-2) / 

J — i 



F n {ix)F m (ix) 4(-l) 

ax = O r , 



cosh 2 (7rx/2) ' h "n(2n+l) 



by a method using the Fourier transform, which we reproduce for a more general class of 
polynomials. The factor (— l) n on the right hand side of (1.2) does not matter, since it 
can be removed by rescaling F n by a factor i n . This factor is necessary in order to make 
the polynomial F n real- valued for imaginary argument. The orthogonality relation (1.2) 
is also derived by Hardy [18, §8] as an example of a general approach to some orthogonal 
polynomials using the Mellin transform, which is equivalent to Bateman's proof of (1.2); 
see also remark 3.2(h). 

The Bateman polynomial F n has been generalised by Pasternack [26] in 1939. He defines 
the polynomial F™ by 

(1.3) F?( 4-) cosh"™- 1 a; = cosh"" 1 " 1 xPJtanhx), 

\ax/ 

for m E C\{— 1}, which reduces to Bateman's polynomial F n in case m = 0. Here we 
use ^cosh~ 1_m x = — (m + 1) cosh -1-771 x tanhx. The case m = 1 already occurs in 
Bateman's paper [10, §4], see also [12]. 

For the polynomials F™ Pasternack derives explicit expressions in terms of hyperge- 
ometric series, generating functions, three-term recurrence relation, difference equations 
and integral representations much along the same lines as in [9], but he does not prove 
orthogonality relations for F™. In particular, Pasternack proves, cf. [26, (10.2), (10.5)], 

r°° l 4 r 1 

/ r(-(m+l + z)) F™(iz)F™(-iz)dz = C (1 - x 2 ) m P n (x)P p (x) dx 

J-oc * J-1 

for some explicit constant C. So if Pasternack would have replaced the Legendre polyno- 
mial P n in (1.3) by the Gegenbauer polynomial P^ n,m \ i.e. the polynomials orthogonal 
on [—1, 1] with respect to the weight function (1 — x 2 ) m , cf. (2.2), he would have obtained 
a one-parameter subclass of the two-parameter continuous symmetric Hahn polynomials 
introduced by Askey and Wilson [4] in 1982. 

Bateman obtains in [12, (3.3)] the orthogonality relations 

f°° F™(ix)F- m (ix) 2(-l) n mn 

(1.4) / — r — -dx = o np — , — 1 < m < 1. 

J -oo cosh(7rx) + cos(m7r) ' 7r(2n + 1) sin7rm 

for Pasternack's polynomials, which reduce to (1.2) for m = 0. The right hand side 
of [12, (3.3)] is not correct. Although explicit expressions for F™, and hence for the 
leading coefficients of these polynomials, were known to Bateman at that time he does 
not rewrite the orthogonality relation (1.4) as the orthogonality relation (1.5) for the 
Pasternack polynomials. 
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In 1956 Touchard [33, §§13, 14] considers orthogonal polynomials associated with the 
Bernoulli numbers. He derives a three-term recurrence relation, orthogonality relations 
and determines the value at — |. In the following paper in the Canadian Journal of 
Mathematics Wyman and Moser [36] give an explicit expression for these polynomials in 
terms of a hypergeometric 4i 7 3-series. A year later Brafman [13] gives another expression for 
these polynomials and derives generating functions for these polynomials. But it is Carlitz 
[14] who notes that the polynomials introduced by Touchard are the same as Bateman's 
polynomials F n defined by (1.1). In 1959 Carlitz [15] replaces the Bernoulli numbers by 
certain numbers involving the Bernoulli polynomials at a point A, where the case A = 
corresponds to the Bernoulli numbers, see also Chihara [16, Ch. VI, §8]. Carlitz shows that 
the corresponding orthogonal polynomials are Pasternack's polynomials (1.3), for which 
he gives the orthogonality relations, see also [4] , 

F™(ix)F™(ix) (-l)»2(l-m) n mn 

r — -dx = o np - - , — 1 < m < 1, 

cos(7rm) + cosh(7nr) ' 2n + 1 n (1 + m) n sin7rm 

where we use the notation of (2.4). The same remark as for (1.2) applies here. The case 
m = of (1.5) is (1.2). From proposition 3.1 we see that (1.5) also yields orthogonal 
polynomials, after a suitable renormalisation, for m G z'R. From either (1.4) or from the 
fact that the weight function in (1.5) is even in m we see that F™ is a multiple of F~ m . 
By comparing leading coefficients, cf. §2, we obtain (1 + m) n F™{x) = (1 — m) n F~ m (x), 
see §2 for the notation. The case m = \ of (1.5) has been given by Hardy [18, §§4-7] in 
1940. 

Carlitz explicitly calculates the moments corresponding to the orthogonality measure 
(1.5) in terms of the Bernoulli polynomials, cf. [15, §6]. This result has already been 
obtained by Stieltjes [29, §5] in 1890 by developing ip(x + b) — ip(x + 1 — 6), ip(x) denoting 
the logarithmic derivative of the gamma function T(x), in powers of x~ x and in terms of 
a continued fraction. This means that Stieltjes gives the three-term recurrence relation 
for the orthogonal polynomials for which the moments of the orthogonality measure are 
given in terms of the Bernoulli polynomials. Stieltjes includes this example in his famous 
memoir "Recherches sur les fraction continues" , in which he also gives the corresponding 
integral representation, cf. [30, §86]. 

In 1982 Askey and Wilson [4] introduce orthogonal polynomials, which generalise the 
orthogonal polynomials introduced by Bateman, Pasternack, Touchard, Hardy and Carlitz. 
These polynomials are orthogonal on R with respect to the measure \T(a + ix)T{^ + ix)\ 2 
(a, 7 > or a = 7 and K(a) > 0) and are nowadays known as the symmetric continuous 
Hahn polynomials. Atakishiyev and Suslov [6, §3] have shown that this is not the end 
of the story and have introduced the continuous Hahn polynomials which have one extra 
parameter, see also Askey [2]. These polynomials are orthogonal with respect to the weight 
function \T(a + ix)T(-y + ix)\ 2 with 9fc(a), 3%) > and 3f(a) = -$(7). 

The goal of this paper is to show that Bateman's approach can be used to prove the 
orthogonality relations for the continuous Hahn polynomials by only using the Jacobi 
polynomials and the Fourier transform in a similar way as Bateman [11] did in order to 
prove (1.2). The orthogonality relations for the continuous Hahn polynomials are not 
new, but this paper shows that Bateman, Pasternack and Hardy could have found these 
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orthogonality relations if they had pursued their approach somewhat further. Moreover, 
this point of view on the relationship between Jacobi polynomials and continuous Hahn 
polynomials gives an intrinsic explanation for the occurrence of the Jacobi polynomials 
in Atakishiyev and Suslov's proof [6, §3] of the orthogonality relation for the continuous 
Hahn polynomials. 

We do not exactly know why Bateman and Hardy did not proceed to find the continuous 
Hahn polynomials as early as the 1930's and 1940's. The following explanation has been 
communicated to me by Richard Askey, whom I thank for letting me reproduce his view 
on this matter here. As to Hardy, we know that he kept the special functions, when 
needed, as simple as possible and that he only used special functions when he had to. So 
going beyond the Legendre polynomial was no option to Hardy. Bateman, as most of his 
contemporaries, thought of hypergeometric series as a function of the power series variable 
z, cf. the definition in §2. But the argument of the Bateman-Pasternack polynomial occurs 
in one of the parameters of a hypergeometric series, cf. remark 2.2(i), and not in the power 
series variable as is the case for e.g. the Jacobi polynomials, and thus they were not in the 
line of thought at that time. This also explains why Rice [28] looked for an appropriate 
generalisation of the Bateman and Pasternack polynomial by introducing a variable at the 
power series spot. However, Rice did not obtain orthogonal polynomials in this way. 

There are more orthogonal systems involving orthogonal polynomials that are mapped 
onto each other by the Fourier transform, or by another integral transform such as the 
Mellin and Hankel transform. The best known example of this are the Hermite functions, 
i.e. Hermite polynomials multiplied by e~ x / 2 , which are eigenfunctions of the Fourier 
transform. For more examples we refer to Koornwinder [22], [23] and to the integral 
transforms of the Bateman project [17]. It is however important to note that in the 
derivation presented here we do not use orthogonal systems, but biorthogonal systems 
involving Jacobi polynomials. This gives the possibility to introduce the necessary extra 
degrees of freedom. Moreover, the result here is not a special case of Koornwinder's result 
in which the Jacobi polynomials are mapped onto the Wilson polynomials by use of the 
Jacobi function transform, cf. [22], [23]. 

A striking aspect of (1.1) and (1.3) is that the argument of the orthogonal polynomial 
is a differential operator. Badertscher and Koornwinder [7], see also [23], have given 
group theoretic interpretations for several identities involving orthogonal polynomials of 
differential operator argument. In these cases these differential operators are acting on 
spherical functions on Riemannian symmetric spaces, which are usually more complicated 
special functions than just cosh -1 x. A related paper in this direction is [21]. 

The organisation of this paper is as follows. In §2 we derive the Fourier transform 
of certain Jacobi polynomials in terms of continuous Hahn polynomials and we discuss 
some applications. We also give the extension of (1.1) and (1.3), and we show how some 
properties of the continuous Hahn polynomials can be derived from properties of the Jacobi 
polynomials. In §3 we prove the orthogonality relation for the continuous Hahn polynomials 
from Parseval's identity for the Fourier transform. 

Acknowledgement. The work for this paper was initiated by the referee report for [21], in 
which the papers by Bateman and Pasternack were mentioned. I thank this referee for 
drawing my attention to these papers. 
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2. Fourier transform on Jacobi polynomials 
The gamma function has been introduced by Euler in 1729 and is defined by 



poo 

\z) = / t z - x e- l dt, U(z) > 0. 
Jo 



The fundamental recurrence relation F(z + 1) = zT(z) follows by integration by parts. A 
closely related integral is the beta integral; 

(2.1) B(a, (3) = f t a -\l - tf- 1 dt = l^^l , »(a), > 0. 

Jo T(a + (J) 

The first proof of this result is given by Euler in 1772. More information on proofs for the 
beta integral and related integrals and sums as well as references to the literature can be 
found in the papers by Askey [3] and Rahman and Suslov [27]. 

The Jacobi polynomials Pn a,/3 \x) of degree n in x are the orthogonal polynomials with 
respect to the beta measure shifted to the interval [—1, 1]; 

(2.2) J (l-x) a (l + x) f3 P^ l3 \x)P^' l3 \x)dx = 

r(n + a + l)T(n + f3 + 1) 
n ' m 2n + a + [3 + l n\ T(n + a + [5 + 1) ' 

for 3?(a) > —1, 9R(/3) > —1. The orthogonality relations are ususally stated for a > — 1, 
(3 > — 1 , but they remain valid under this more general condition on the parameters a and 
p. The weight function is positive if and only if a and (3 are real. An explicit expression 
for the Jacobi polynomial Pn a, ^\x) is given by a terminating hypergeometric series; 

(2.3) p(«.«(x) = <^ s f 1 + " + 1 ; i^) , 

where the terminating hypergeometric series is defined by 

/ -n,ai,... ,Op \ _ A {-n) k {a 1 ) k ...{a p ) k z k 



(2.4) (o) fc = o(o + l)(o + 2) . . . (a + k - 1) = ^""V^ , k E Z+. 

1 (aj 

More information on Jacobi polynomials can be found in Szego's book [31, Ch. 4]. 

Let us calculate a Fourier transform involving Jacobi polynomials. Rewrite the Fourier 
transform 



(2.5) 



/oo 
e~ lxz {l - tanhx) Q (l + tanha^P^'^tanha;) dx 
-oo 

= J 1 (1 - t) a ~ 1+ ^ z (l + tf^-^p^it) dt 

= 2 a+/3-l f u *-^H™(l- U )P-l-¥zp("lf)(l-2u)du 

Jo 
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using the substitutions t = tanhx, t = 1 — 2u. Use that 

4- = (cosha;) -2 = (1 — tanhx)(l + tanhx), - — j^jj^ _ e ~ 2x _ 
dx K ' K A ' 1 + tanhx 

Note that x i-> (l-tanhx) a (l+tanhx)^ 7,5) (tanh;r) G L 2 (R)nL 1 (R) for 3?(a), 3?(/3) > 0. 

In (2.5) we use the explicit series representation (2.3) for the Jacobi polynomial and the 
beta integral (2.1) to see that (2.5) equals 

t ( - n) t ( ;i^ +1)t ^ h 



fc=0 



k 



(2.6) = y^-l (7 + l)n + j fa g _ ^ (-")*(" + 7 + *+ + * fa >> 

= 2»W- <l±^B(a + i«, /J - i«) 3 F 2 f " + 7 +/ + 1 ' ° + , A . 

n\ 2 2 V 7 + l,Q! + /5 / 

The identity obtained in this way can also be found in Erdelyi et al. [17, 16.4(3)]. 
The 3F2-series in (2.6) is a continuous Hahn polynomial defined by, cf. [2], 

. n (a + c) n (a + d) n f-n,n + a + b + c + d-l,a + ix \ 

(2.7) p n (x;a,b,c,d) =i 3 F 2 ;1 . 

nl \ a + c, a + d J 

So we have proved the following lemma. 

Lemma 2.1. For z£R, 3R(a), 3R(/3) > and -7 p toe /icwe 

/oo 
e" lX2 (l - tanh;r) Q (l + tanha^P^'^tanhz) dx 
-00 

1 (a + p + n) 

where Pn ,S \x) is a Jacobi polynomial defined by (2.3) and p n (x; a, 6, c, d) is a continuous 
Hahn polynomial defined by (2.7). 

Remark 2.2. (i) An equivalent way of proving the lemma is by establishing 

(2.8) p n (---^-; a, 8 - p + 1, 7 - a + 1, p){(l - tanhx) a (l + ttrnhxf} = 

2 dx 

i n (a + /3) n {(l - tanhaO a (l + tank xf} P<?' s) (tank x) 

and applying the Fourier transform T to it and using T(f')(z) = izJ r (f)(z). Equation 
(2.8) can be proved from 
(2.9) 
d 

— [(1 -tanhx) a (l + tanhx) /3 ] = (a + [5 + (a - 0) tanhx) [(1 - tanhx) a (l + tanhx)' 3 ], 

(July 

(a + --^-) r [(l - tanhx) a (l + tanha^l = 
2 dx 

2" r (l -tanhx) r (a + /3) r [(l -tanhx) a (l +tanhx) /3 ], 
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which in turn can be proved by induction with respect to r G Z + . Special cases of (2.8) 
are (1.1) for a = (3 = |, 7 = 8 = and (1.3) for a = f3 = \{m + 1), 7 = 5 = for the 
Bateman polynomial and its generalisation by Pasternack. So we have 



( -n,n+l,±(l + m + x) \ 

" 3 2 V 1, m + 1 ' ; 



for the Bateman (m = 0) and Pasternack polynomials defined in (1.1) and (1.3). The case 
a = f3 = \ + I A, 7 = 5 = — \ of (2.8) has been observed by Koornwinder [23, §2]. 

(ii) Applying the inverse Fourier transform and taking n = gives the Fourier transform 
of z 1— > F(a + iz)T((3 — iz), which is closely related to the orthogonality measure for the 
Meixner-Pollaczek polynomials, cf. Askey [3]. 

Lemma 2.1 can be used to obtain identities for the continuous Hahn polynomials from 
identities satisfied by the Jacobi polynomials. As a first example we start with the following 
generating functions for the Jacobi polynomials, cf. [20, (1.8.7), (1.8.6)], 

(1 _ tn - 4 - ljFl (l(7 + * + l),l(7 + * + 2) ; ^J = 

n=0 17 + L>n 

and 

( - [x-\)t \ ( - (x + \)t \ ^ Pt?- S \x)t" 

^Ki + x'—T- ) " Fl {s + r^-)-2-. h + iU6 + i) n 

A straightforward manipulation using lemma 2.1 proves the following generating func- 
tions for the continuous Hahn polynomials. The first of these generating functions is also 
contained in [20, (1.4.8)]. We get 

n ,,_ a _ /3 _ 7 _ 5 _ 1 f±(a + + 7 + 6-l),±(a + + 7 + 6),a + iz, -At \ 

™ ( a + P + 1 + 5-l) n (z . a6 „ 0) 
^ (a + j 9) n (a + 7 )„ K , ) ^^' a ' () ' 7 '^' 

which is the generating function used by Bateman [9, §3] and Pasternack [26, (2.2)], and 

(t/i) n Pn (z;a,6,>y,P) (-t)H k (a + iz) p ((3 - iz) k 



(7 + <*) n (S + (3) n {a + (3) n ~ ^ ^ p! (7 + a) P fe! (5 + /%(« + ' 

The last series can be rewritten as a hypergeometric series in two variables with arguments 
x = -t,y = t, cf. Appell and Kampe de Feriet [1, Ch. IX, p. 150, (29)]. 



8 



H.T. KOELINK 



As a second example we derive two relations between three continuous Hahn polyno- 
mials from identities for Jacobi polynomials. First use J r (f'){z) = izJ r {f)(z), (2.9) and 
d_ P {i,5)^ = i( n + 7 + 5 + l)p^+ 1 >* +1 )(a;) J c f. [31, (4.21.7)], and straightforward calcu- 
lations to get after some rewriting 

(a + (3 + n)izp n (z; a, 8, 7,/?) = (a + (3)(a + iz)p n (z; a + 1,5,7- 

+ z(n + a + /3 + 7 + <5 - l)(a + iz)((3 - iz)p n -\{z\ a + 1, 8, 7, (3 + 1). 

Another classical identity for the Jacobi polynomials, cf. [31, (4.5.4)], 

(n + 7 + 1)^ 7 ' 5) (x) - (n + l)P^f (*) = ^ (2n + 7 + 5 + 2) (1 - x)P^ +1 ' <5) (x) 
leads to 

(2n + a + /3 + 7 + <5)(a + zz)p n (z; a + 1, 5, 7, /?) = 

(a + /3 + n)(n + 7 + a)p n (z;a, 5,7,/3) + i(n + l)p n+1 (z; a, 8, 7, (3). 

It is also possible to use two or more identities for the Jacobi polynomials in order to 
obtain identities for continuous Hahn polynomials. As an example we indicate how the 
three-term recurrence relation for the continuous Hahn polynomials can be derived, cf. 
Pasternack [26, §5]. Let p n , cf. (2.8), be defined by 

(2.10) p n (-^-){(l-tanhx) a (l+tanha:) /3 } = {(1 -tanhx) Q (l + tanha:) /3 }P 7 ( 7 ' 5) (tanhx). 
Differentiate this identity once more to get 

-l-Pn(-^-){(l - tanha:) a (l + t&nhx) 13 } = {(1 - tanha:) a (l + tanhxf} 
dx ax 

x [(a + f3 + (a - f3) tanhx)P 7 ^ 7 ' ,5) (tanha:) + (1 - tanh 2 x) -|- (tanhs)] . 
In the term in square brackets we use 

(i - * 2 )^|r^) = AnP^ix) + B n p<?>*>(z) + c n p™(x), 

cf. Szego [31, (4.5.5)] for the explicit values of the constants, and the three-term recurrence 
relation for the Jacobi polynomials, cf. [31, (4.5.1)], to get only Jacobi polynomials of 
degree n + 1, n and n — 1. Recalling (2.10) we find the three-term recurrence relation for 
the continuous Hahn polynomials p n . For the explicit values of the coefficients we refer to 
[20, (1.4.3)]. 
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3. Orthogonality for the continuous Hahn polynomials 
The set of functions x \— > (1 — tanhx) a (l + tanh x) 13 pji 2a+1,2f3+1 ^ is an orthogonal basis 



of L 2 (M) and by lemma 2.1 it is mapped by the Fourier transform onto the set of functions 
z i— > T(a + ^iz)T((3 — ^iz)p n (^z; a, P,a, P). Since the Fourier transform is isometric we 
obtain the orthogonality relations for the continuous symmetric Hahn polynomials, cf. [4] , 
see also [22], [23], but we can do better as follows. 

From the Parseval identity 2tc J r f(x)g(x)dx = J R (^Ff) (z) (Fg) (z)dz for the the Fourier 
transform for /, g G L 2 (R) we obtain 

/oo 
(1 - tanhx) a+a (l + tanhx) /3+6 Pi 7 ' ,5) (tanha:)Pi c ' d) (tanhx) dx 
-oo 

(3.1) = im -n 2a+a+ p +b -2 [°° r(a + jiz)T(J) - \iz) Tja - liz)T(b + \iz) 



x p n (z/2; a, 8 — (3 + 1, 7 — a + 1, P)p m (z/2; a 7 d — b + 1, c — a + 1,6) dz 

for dt(a), 9?(/3), 3?(a), 3?(6) > 0. Next we restrict the parameters in the left hand side of 
(3.1) such that we get the orthogonality relations (2.2) for the Jacobi polynomials. So 
we take 3?(a + a) > 0, + 6) > 0, which is already satisfied, and r y = c= a + a — 1, 
5 = <i = /3 + 6 — 1. For these choices we put again t = tanhx to see that the left hand side 
(3.1) equals zero for n 7^ m. The square norm follows from (2.2) and so we get 



(3.2) 



/oo 
T(a + iz)T(/3 — iz)T(a — iz)T(b + iz)p n (z; a, b, a, f3)p m (z; a, j3, a, b) dz 
-00 

- 2n5 r< " a + 13 + n } T ( a + b + n ) r ( n + " + a ) r (^ + P + b) 
71 n,m n\(2n + a + P + a + b- l)T(n + a + P + a + b-l) 



subject to the conditions K(a), 9£(/3), 3?(a), 3?(6) > 0, since 7 = c = a + a — 1, 8 = d = 
P+b—1. If we replace p m in (3.2) by lc{p rn )z rn then (3.2) remains valid for < m < n. Here 
lc(p m ) denotes the leading coefficient of p m (z;d,P,a,b), and since this equals the leading 
coefficient of p n (z; a, 6, a, P), cf. (2.7), we may replace p m (z; a, j3, a, 6) by p m (z; a, 6, a, /?). 
So we have proved the orthogonality relations for the continuous Hahn polynomials stated 
in the next proposition from the Parseval identity for the Fourier transform and from the 
explicit knowledge of the Fourier transform of the Jacobi polynomial described in lemma 
2.1. 

Proposition 3.1. The continuous Hahn polynomials defined in (2.7) satisfy 
1 f°° 

(3.3) — / r(ct + iz)T(P — iz)T{a — iz)T{b + iz)p n (z; a, 6, a, P)p m (z; a, 6, a, P) dz 

r(a + /3 + n)r(a + b + n)T(n + a + a)T(n + P + b) 
n,m nl(2n + a + P + a + b-l)F(n + a + P + a + b-l) 

for K(a), 3?(a), 3?(6) > 0. 

Remark 3.2. (i) The weight function is positive for a = a, 6 = P, or for a = j3, a = 6, 
which follows from the invariance of (3.3) under interchanging a and b or P and a. 
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(ii) The case n = m = of (3.3) is Barnes' first lemma from 1908, see e.g. Bailey [8, 
1.7] or Whittaker and Watson [34, §14.52] for the original proof by Barnes. An equivalent 
proof of Barnes' first lemma as given here can be found in Titchmarsh's book [32, (7.8.3)], 
where the Mellin transform is used instead of the Fourier transform. Proposition 3.1 can 
be obtained in a similar way using the Mellin transform if we use the Jacobi polynomials 
of argument (1 — x)/(l + x). We obtain an orthogonal system on [0, oo). So we start with 
the Mellin transform 



and the Parseval formula for the Mellin transform gives proposition 3.1. See Hardy [18, 
§8] for this approach to Bateman's polynomial F n . In [24, prop. 3.1] Koornwinder shows 
that the Laguerre polynomials are mapped onto the Meixner-Pollaczek polynomials by the 
Mellin transform. The Mellin transform of the underlying measures is given in Titchmarsh 
[32, (7.8.1)]. This is a limiting case of (3.4). Replace in (3.4) x by x/8 and (3 by 5^ with 
9?(£) > 0, 1^(01 < 71 an d let 5 — > oo. Similarly we can obtain the analogue of lemma 2.1 
with the Laguerre and Meixner-Pollaczek polynomials by a suitable limit transition. 

(iii) The result (3.3) in this form has been proved first by Askey [2] using Barnes' first 
lemma and the Chu-Vandermonde summation formula for a terminating 2i 7 i-series of unit 
argument. Before that Atakishiyev and Suslov [6], see also [25, §3.10.3.2], proved (3.3) 
in case of a positive weight function. The method employed by Atakishiyev and Suslov 
uses Barnes' first lemma, which is rewritten in terms of the beta integral so that the 
orthogonality relations of the Jacobi polynomials can be used. 

(iv) Another proof of proposition 3.1 using symmetry in the parameters a, 6, a, (5 is 
given by Kalnins and Miller [19, §3]. They also give a proof of Barnes' first lemma in this 
way. 

(v) The continuous Hahn polynomials are not on the top shelf of the Askey-scheme of 
hypergeometric polynomials. The most general hypergeometric orthogonal polynomials 
with a continuous weight function are the Wilson polynomials, cf. Wilson [35], which 
have four degrees of freedom. The continuous Hahn polynomials can be obtained from the 
Wilson polynomials by a suitable limit process, cf. e.g. [20, §2.2]. 

The orthogonality relations stated in §1 are special cases of (3.3). In particular, the 
orthogonality (1.5) for the Pasternack polynomials follows by taking a = f3 = \(1 +m), 
a = b = |(1 — m). This shows also that we have a positive weight function for — 1 < m < 1 
or m G z'R, which are equivalent for m and — m. To see this we have to use the reflection 
formula T(z)T{l—z) = n sin _1 (7rz) and some straightforward manipulations on goniometric 
and hyperbolic functions. It should also be noted that taking the same values for the 
parameters in (3.2) gives Bateman's (bi)orthogonality relations (1.4). 



(3.4) 




T(a-iX)T(P-iX) 
T(a + (3 + n) 



i n p n (-X;a,S - /3+ 1,7 - a + 1,/?) 
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